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Abstract 


T lepton is one of three chareged leptons in nature, the measurements of its mass 
have been performed since its discovery. The present relative accuracy is already 
at the level of 1074; more factors are still being studied in order to increase the 
accuracy. However, the available techniques for analysis and expectable luminos- 
ity from ete collider indicate that the precision upper limit of 7 mass is almost 
reached, which means that brand new approaches should be looked for if the great 


improvement is yearned for. 
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1 Introduction 


The history of 7 mass (m,) measurement has forty years. In the first experimental 
paper of 7 lepton [1], m, is estimated to have a mass in the range from 1.6 to 2.0 GeV. 
Since then many experiments have been performed to measure m. [2, 3, 4, 5, 6, 7, 8, 9, 
10, 11, 12, 13, 14, 15, 16, 17], the results are displayed in Fig. 1. 

The measurement results of m, in the 21 century are summarized in Table 1, where two 


results were acquired using the method of pseudo-mass while the others using the method 
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of threshold scan. For the pseudo-mass method, the huge amount of data from B-factories 
are employed [14], good statistical accuracy is acquired, but large systematic uncertainty 
exists, which is mainly due to the absolute calibration of the particle momentum. For the 
threshold-scan method, the value of m, was extracted from the dependence of production 
cross section on beam energy. In the KEDR experiment [15], both resonant depolarization 
technique and Compton backscattering technique [22] are used to determine the beam 


energy. All these techniques greatly decrease the uncertainty of beam energy. 


Table 1: Measurement results of m, in the 21 century. 


Measured m; Year | Exp. Group | Data sample | Method 
(value+statistic+systematic) 

1776.91 + 0.127038 2014 | BESIII [17] | 23.26 pb^! | Threshold-scan 
1776.68 + 0.12 + 0.41 2009 | Babar [16] | 423 fb! Pseudo-mass 
1TOO RI OS +0.15 2007 | KEDR [15] | 6.7 pb! Threshold-scan 
1776.61 + 0.13 + 0.35 2007 | Belle [14] 414 fb! Pseudo-mass 


Although the accuracy of results from above two methods are at the comparable level, 
it is obvious that the systematic errors already dominate for the pseudo-mass method, the 
enhancement of luminosity is already unappealing. For the threshold-scan method, both 
the statistic and systematic errors still seem to have room for the further improvement. 


Therefore, the subsequence discussions will be devoted to the threshold-scan method. 


2 Statistic aspect 


By the virtue of optimization theory [23], the relation between the absolute error of m. 


(denoted by u) and the total luminosity (denoted by L) is given by the following formula 
= 
vL’ 


where A is a constant depending on the cross section, the derivative of cross section to 


(1) 


m,, and some other quantities [23]. 
An effective error variation that is called L-profit, is introduced and defined as 
Au du A 3 
= li ———/|=-—_=—[2. 2 
T= ADS | A dL 2 (2) 
The minus in the above definition indicates that when L increases u always decreases. 
This quantity discloses the variation of u for a unit luminosity. Although it is reasonable 


to expect the reduction of error when the luminosity accumulates, but L-profit will also 


be small for larger L. Refer to Fig. 2, when u is considerably large, the enhancement of 
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Figure 1: m, measured in the last and this century. In (a) the red line indicates the 


(b) m. measured in the 21 century 


average value of m, in PDG2000 [18]: m, = 1777.03 


measured value from BES in 1996 is also plotted in (b) but with blank blue circle for 
distinction. The average value of m, in PDG2015 [19]: m, = 1776.86 + 0.12 MeV is also 
It should be noted that since PDG1996 [20, 21], the 


results from experiements performed before 1990 were removed except for the result of 


indicated by the red line in (b). 


DELCO. 


MeV. For comparison, the 


L can effectively lead to the prominent improvement of u. But step by step, when u is 
already accurate enough, a great mount of L has to be consumed for reducing the error 
while the improvement is very limited, which is rather uneconomical. 

The above analysis is not the whole story. As a matter of fact, the increment of L 
means that more machine time must be used, more data must be dealt with, and more 
analysis must be performed. Moreover, the improvement of u means that a lot of work 
concerning systematic uncertainty analysis must be finished, which in turn leads to the 
necessary demand for improvement of Monte Carlo simulation. A great amount of details 
has to be taken into account, sometime the upgrade of software is indispensable as well. 
Intuitively the work consumed for error improvement increases exponentially, and this 
is even more real when error is already small enough!. Under such a consideration, a 
work-factor D is designated as 
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B 
D =e% = evt. where C 2 4. (3) 


Here the constant B will be related to concrete analysis procedure. Similar to Eq. (2), 


l'There is no induction model or deduction proof for this exponential assumption, since the work 
done for a certain target is usually hard to be quantified. But for such an assumption, an interesting 
fact is worthy of mentioning which may be edificatory for understanding the exponential accretion. In 
the standard model the contribution to the anomalous magnetic moment ae comes from three types of 


interactions, electromagnetic, hadronic, and electroweak: 
ae = a, (QED) + ac(hadronic) + ae(elctroweak). 


The QED contribution can be evaluated by the perturbative expansion in a/7: 


ae(QED) = Y: (S) a, 


n=1 


where aC"? is finite due to the renormalizability of QED and may be written in general as 


af = AP? + AP” (me/my) + AY” (mu/mr) + AS” (me/Mp, my mz) 


to show the mass-dependence explicitly. Agr) can be calculated by QED theory perturbatively with the 


increase of the order n. It is noticeable that the number of Feynman diagrams increase even more rapidly 
than the exponential increase with respect to n. The following table shows the interesting comparison [24, 
25]. 


Number of Feynman diagrams (N.F.D.) for Av") 


order n=1 n=2 n=3 n=4 n=5 
2n 2 4 6 8 10 
AQ9 A0 A40 49 49 409 
N.F.D. 1 5 72 891 12672 
eYn-1) 1 7.39 54.60 403.43 2980.96 


another effective error variation that is called D-profit, is introduced and defined as 


Au du/dL A? 
= lim |- -- a2 ee. 4 

£ Jim. d oe UE Y 4) 
If refer to Fig. 2, it is fairly clear that the decrease of £ with L is violent, which indicates 
the increase of luminosity will play smaller and smaller effect on the improvement of m». 
In another word the error improvement of m, in e*e- collider is no room using current 
method. 
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Figure 2: The relation between luminosity and error. The black solid line denotes the 
variation of u against L; the red dash line the variation of v against L; the blue dotted 


line the variation of € against L. In calculation A = 1 and B = 1 are adopted. 


To have a concrete impression of above analysis, we compare two pedagogical experi- 
ments, the first one intends to decrease the error from 2 MeV to 1 MeV, while the second 
one intends to decrease the error from 0.2 MeV to 0.1 MeV. The coefficients A and P are 


chosen to be 1 for the unit of MeV. The ratio of aforementioned quantities is defined as 
R, = (5) 


where superscript i indicates the first (i = 1) or the second (i = 2) experiment; subscript 
a means the value after experiment, while b the value before experiment. The symbol 
X denotes the quantity L, v, D, and € respectively. The evaluations are presented in 
Table 2. 

Based on Table 2 we know that if the error is improved to be one half of the previous 


experiment, four fold luminosity is required, and luminosity profit is one eighth. Actually 
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Table 2: Evaluations for pedagogical experiments. 


Ratio value First Exp. Second Exp. 
2— 1 MeV 0.2 0.1 MeV 
Hye (uu 4 4 
R, = (ua ug)? 1/8 1/8 
Rp = eUfus-lfu) 1.65 148 
Re = (us/ up)? - e 71/uo) 0.15 1.69 x 107? 


according to quantities L and v, as long as the improvement proportion is the same, Rz 
and R, are the same. This also implies that these two quantities can reflect neither the 
realistic working strength nor working profit. Turning to another two quantities Rp and 
Re, the situation becomes rather different. Although the improvement proportion is the 
same for the two experiments, the work have to be done for the second experiment is 
ninety times more than that for the first one, while D-profit for the second one is merely 
one ninetieth of the first one! This reminds us of the statement, that is, when u is already 
accurate enough the further increase of luminosity is uneconomical and ineffective. 

A remark is in order here. During the study of m, scan, BESIII collaboration have 
make methodical and meticulous studies on optimization process, to guarantee the ef- 
fective usage of available luminosity [23, 26, 27, 28, 29]. Nevertheless, this effort is not 
helped with the low Re issue. 


3 Limit due to systematic errors 


For BESIII collaboration, before the actual scan of m,, many studies have been made, 
one of which is the systematic estimation [26] as listed in Table 3. It can be seen that 
the error due to energy scale dominates. To decrease such an uncertainty, starting from 
year 2007, a high accuracy beam energy measurement system (BEMS) located at the 
north crossing point (NCP) of BEPC-II was designed, constructed, and put into the 
commissioning at the end of 2010 [30, 31, 32, 33]. The launching of the system is excellently 
well, two days are utilized to perform v" resonance scan. The mass difference between 
the PDG2010 value [34] and the measured result by BEMS is 1 + 36 keV, the deviation 
of which indicates that the relative accuracy of BEMS is at the level of 2 x10? [32]. 

In the actual data taking process, it is found that the measurement accuracy of BEMS 
is sensitive to the running status of accelerator [35]. At the end of 2011, a test scan of m; 
was performed. The integrated luminosity around of 23.26 pb^! is collected in the vicinity 
of r-pair threshold, together with 1.5 pb ^! and 7.5 pb ^! data sets for the J/1» and w’ 


resonance scans respectively, by virtue of which the systematic uncertainty due to beam 


Table 3: Systematic uncertainty estimation for m, measurement from Ref. [26]. 


Source om, dm,/m, 

(107? MeV) (1079) 
Luminosity 14.0 7.9 
Efficiency 14.0 7.9 
Branching Fraction 3.5 2.0 
Background 1.7 1.0 
Energy spread 3.0 1.7 
Theoretical accuracy 3.0 1.7 
Energy scale 100 56.3 
Summation 102 57.5 


energy is determined. The mass of 7 lepton is measured from a maximum likelihood fit 
to the 7-pair production cross-section to be m, = (1776.91 + 0.12*015) MeV [17], among 
which the systematic uncertainty due to energy scale also dominates all relevant terms, as 
shown in Table 4. This tough situation is the same for all other ete~ colliders that adopt 
the threshold-scan method, which can not be circumvented either for the pseudo-mass 
method. 

The crucial issue here lies in the uncertainty of energy scale. For BESIII at BEPCII, 
BEMS was established to control this uncertainty. The working scheme of this system can 
be described briefly as follows [36]: firstly, the laser source provides the laser beam and 
the optics system focuses the laser beam and guides it to make head-on collisions with 
the electron (or positron) beam in the vacuum pipe, here the Compton backscattering 
process happens; after that the backscattering high energy photon will be detected by the 
HPGe detector, the detection capacity of which pins down the accuracy of beam energy. 
For the time being, the limit of relative accuracy for HPGe detector is at the level of 1075, 
which in turn restricts the limit for determination of beam energy?. Another juxtaposing 
restriction at the level of 1075 comes from the long term stability of electric power supply 
of BEPCII [41]. This stability directly controls that of magnets, which determines the 
energy of beam. 


For the present BEMS, the detector calibration is the main factor that restricts the detection accuracy. 
Generally speaking, there are two methods to calibrate the HPGe detector, one is the radiative source 
calibration. Several radiative sources are put nearby the detector and used to provide the simultaneous 
calibration lines during the data taking. The demerit of this method lies in that sometime not suitable 
sources can be found in the wide region of measurement. The other calibration method is to utilize the 
radiative sources from nuclear reactor, which can provide a great amount of calibration lines. But the 
shortcoming of this method consists in that calibration can be not done simultaneously with the data 
taking. However, in principle the calibration of detector can be improved further. 


Table 4: Summary of the 7 mass systematic errors. The meanings of each term can be 
found in Ref. [17]. 


Source Am, (MeV/c?) 
Theoretical accuracy 0.010 
Energy scale re 
Energy spread 0.016 
Luminosity 0.006 
Cut on number of good photons 0.002 
Cuts on PTEM and acoplanarity angle 0.05 
mis-ID efficiency 0.048 
Background shape 0.04 
Fitted efficiency parameter rry 
Total SET 


4 Way out for accuracy improvement 


The preceding analyses obviously indicate that the efforts spent on e*e^ collider are 
not to bring the promising improvement on the accuracy of m+. Therefore, the new ap- 
proach indeed needs to be figured out in order to increase the accuracy of m, significantly. 
Without any plausible precedent or theoretical implication, we would like to look at some 
relevant measurements with high accuracy to fish out some clues for m, experiments in 


the future. 


4.1 Experimental aspect 


According to PDG2012 [42] 


Me 0.510998910 + 0.000000013 MeV — (óm, m, ~ 2.554 x 1075) , 


105.658367 + 0.000004 MeV | (6m,,/m, œ 3.786 x 10-8) , (6) 


My 


Such a highly precise determination of an electron’s mass comes from measuring the 
ratio of the mass to that of a nucleus, so that the result is obtained in m,, (atomic mass 
units). For example, by comparing cyclotron frequencies of electrons and single C9* 
ions alternately confined to the same uniform magnetic field in a Penning trap |43], the 


electron's mass is determined as 


Me = 0.0005485799111(12)m, . 


Then a convert factor 941.494013 + 0.000037 MeV /m,, is used |42, 44] to transform unit 
Ma into MeV. 

As far as the muon's mass is concerned, it is obtained from the muon-electron mass 
ratio as determined from the measurement of Zeeman transition frequencies in muonium 
(u*e- atom). Unfortunately, the life time of 7 lepton is too short to form such an atom, 
so the mass of 7 lepton can not be determined accurately in the similar way. 

Anyway, it is fairly enlightened that we should find a relation of m, with another 
sensitive quantity which can be measured accurately, or their ratio can be measured 


accurately. This might be the new direction for the further accurate measurement of m;. 


4.2 "Theoretical aspect 


Similar to the situation of experiment, theoretically speaking, if a certain relation can 
be found for m,, the excellent accuracy may be expected. Here we mention an interesting 


formula of masses of three leptons [45, 46, 47, 48]: 


2 
Me + m, 4 m, = z(y Me + My + mr)” , (7) 
In the light of which it is easy to find 


VT = 2(v/mz + Vig) + /3(m« + My) + LDV Tie Mg (8) 


Using this formula together with the mass values for electron and muon in Eq. (6) the m. 


can be obtained together with the corresponding error as follows 
m, = 1776.968884 + 0.000058 MeV (ém,/m, zz 3.26 x 1078). (9) 


Unfortunately, the formula (7) itself is not obtained from the first principle of particle 


physics, so this kind of efforts can only be treated as an attempt direction in the future’. 


3Among the mysteries of particle physics, one of the most elusive is the mass spectrum that can't 
be predicted from the elementary principle, which is also the most unsatisfactory feature of the present 
Standard Model. The arbitrariness of mass assignment triggers the decades of studies [37, 38, 39, 40], 
among which several interesting results emerge. A very simple formula (i.e. Eq. (7)) between the masses 
of the three charged leptons was stumbled over by Yoshio Koide at the beginning of 1980's when he was 
working on some composite models of quarks and leptons [45, 46, 47, 48]. Further relevant studies are 
performed [49, 50, 51, 52, 53], although the impact is small, some interesting results are worth of mention, 
one of which is Barut's formula for leptons [51], viz. 


3 N 
-1 4 
mN) =m. (1430 nt) 


where me is the mass of electron and o ~% 1/137 the electromagnetic fine structure constant. The masses 
of electron, muon, and tauon correspond to N = 0, 1, and 2, respectively. 


5 Summary 


In this paper, the issue concerned the upper limit of the error of m, measuremen- 
t is discussed from unconventional point of view. The main conclusions acquired are 


summarized below: 


1. Statistically speaking, the increase of luminosity have no effective effect on the accu- 
racy improvement of m,. Moreover, from point view of efficiency and profitability, 


the input-output ratio will considerably decrease as the accuracy of m, increases. 


2. Systematically speaking, the limit at the level of 107? of both calibration of HPGe 
detector and energy stability of accelerator circumscribes the accuracy limit of m,. 
For experiments at ete collider, this is the insurmountable limit for the systematic 


uncertainty. 


3. The absolutely distinctive train of thought from both theoretical and experimental 


aspects is needed to seek a new direction to increase the accuracy of m,. 
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